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ABSTRACT 


Following success in other areas as a remote sensor, the spacebome microwave 
image radars are assuming a notable position in the problem of ship detection for civilian 
and military purposes. This work will discuss the strong and weak points of Synthetic 
Aperture Radar (SAR) when used for ship detection. First, the thesis gives a brief 
description of SAR fundamentals, image processing and the parameters for ship 
detection. Second, the actual techniques, limitations, errors and some models used for 
ship detection are described. Finally, using a well-known and reliable ship detection 
model (Vachon et al. 1997), tested in the Canadian Ocean Monitoring Workstation and in 
some validation field programs, the new generation of spacebome SARs, mainly 
RADARS AT 2, are analyzed for ship detection capabilities. During the analysis 
parameters like wind velocity, wind direction related to the antenna, satellite incident 
angle and Number of Looks are changed to study their influence on ship detection. 
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1. INTRODUCTION 


A. MOTIVATION 

Having approximately 7,480 km (4650 miles) of coastline and 200 nautical miles 
of territorial sea, the total marginal sea area under the supervision of the Brazilian Navy 
totals the impressive amount of approximately 2,773,000 knf. Maritime control and 
surveillance of this huge area requires all types of tools. To fulfill this task in time, with 
reliable accuracy, airborne and spaceborne microwave imaging sensors emerge as a 
valuable instrument. The advantages of this kind of tool include a regular repeating 
observation cycle, images free of cloud interference, high ground resolution (up to 3m in 
RADARS AT-2) and reliability. Synthetic Aperture Radars (SAR) have these essential 
features and offer a good trade-off between cost and benefit. In the future, use of SAR 
will be supported by a commercial remote sensing imaging boom. The new commercial 
SAR systems can supply images for specific requests with competitive cost and a 
potential for integration within modular ship monitoring systems. An increasing number 
of missions are being planned and executed as outlined in Table 1. 

This thesis focuses specifically on three missions that will fly the next generation 
SAR sensors designed to fulfill various applications including ship detection. These 
sensors will be launched between 2001 and 2003 and include the following: the Japanese 
Space Agency Advanced Land Observing Satellite with the Phase Array type L-band 
Synthetic Aperture Radar (ALOS/PALSAR); the European Space Agency Environmental 
Satellite with an Advanced Synthetic Aperture Radar (ENVISAT/ASAR) and the 
Canadian Space Agency Radar Satellite 2 (RADARSAT 2). 
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Table 1. 


Planned SAR satellites missions [Updated from Olsen et al. (2000)] 


Satellite 

(Owner/Operator) 


ENVISAT (USA) 

Launch 

jan 2002 


RADARSAT-2 
(MDA/CSA) 
Launch 2003 


ALOS (NASDA) 
Launch 2003 


TerraSAR 

(DLR/DASA/MMS/BN 
SC) [4] 


COSMO/SKYMED 


Lights AR 
(NASA/Industry) 


Earth Watch 
(ESA/Industry) 


Description 


ENVISAT will be the next 
satellite to be launched by the 
European Space Agency. The 
C-band ASAR instrument will 
be the first dual-polarised 
orbital SAR with a long term 
presence in space [3]. 


This mission is designed to be 
yet another step closer to a 
fully commercial mission. 
Most of the investment is 
carried by the Orbital 
Sciences Subsidiary, Mac 
Donald Dettwiler, with the 
Canadian Space Agency as an 
anchor customer. The satellite 
will carry a C-band SAR that 
supports viewing to either side 
of the spacecraft track, high 
and low resolution modes and 
polarimetric capability. 


Primarily a land observing 
mission, operating at L-band. 
The wide swath mode 
provides an imaging geometry 
that will provide good target 
to clutter ratios. 


TerraSAR is a commercial 
mission planned by a 
consortium of British and 
German industry. Current 
plans indicated a two-satellite 
system with X and L-band 
capability. 


Key radar parameters 


Swath width 


56-105 Km 
56-105 km 
405 km 
405 km 


Incidenc 
e angles 



Image 

ALTPOL 

Wide 

Global 


Standard 
Wide 
Low Inc 
High Inc 
Fine 

ScanSAR Wide 
ScanSARNarrow 
Std QuadPod 
Fine QuadPod 
Triple Fine 
Ultra Fine Wide 
Ultra Fine Narrow 


Image/Dual pol 

modes 

ScanSAR 


X-band: 
10 km 
20 km 
100 km 
L-band: 
40 km 
60 km 
200 km 


Spotlight 

Stripmap 

ScanSAR 


Quadpol 

Stripmap 

ScanSAR 


This IS presently a mission 3 km <45 1 m Spotlight 

concept developed under 30 km < 45® 3 m Stripmap 

funding from the Italian Space 120 km < 45® 25 m ScanSAR 

Agency. The SAR system 

considered is the so-called 

SAR2000, an X-band radar 

capable of stripmap, spotlight 

and ScanSAR imaging. 


This mission has been planned for some time by NASA, where an industry partner 
has been sought. Currently, discussions are being held with oil industry, to determine how a 
mission can be designed to meet that industry sector’s requirements as well as NASA’s science 
goals. The proposed radar is an L-band system with moderate resolution. 


The ESA Earth Watch programme is intended to be ESA’s contribution towards a 
European Earth observation and monitoring system, in cooperation with private industry and 
organizations such as EUMETSAT. Industry proposals for mission concepts have been 
solicited, and several SAR missions are included. Further selections have not been confirmed 
to date. 












































This thesis presents a technical review of one of these spaceborne SARs - the 
Canadian RADARSAR 2. This choice is due to the availability of commercial imagery, 
access to public information on the sensor, the success of the previous mission 
(RADARSAT-1), the set of technical improvements, and the suitability for the Brazilian 
Navy Project. Issues like ship tracking and target recognition that result from ship 
detectability improvements will be mentioned but not pursued in detail in this thesis. 
Expected results with new polarization improvements will also be mentioned but not 
analyzed in detail. This thesis will describe SAR operations and the physics of the ship 
detection problem. 

B. THESIS OUTLINE 

The first chapter describes the motivation and focus of this thesis. Chapter II 
presents some necessary SAR fundamentals and background needed to understand the 
ship detection problem. Constraints, limitations and explanation of system modes 
precedes the SAR imaging process in order to present a complete view of data acquisition 
through final image creation. Chapter III deals ^ecifically with all facets of the ship 
detection problems, limitations, errors and techniques. The fourth chapter is the thesis 
core. All the results of quantitative improvements for the next family of spaceborne 
SARs are included in Chapter IV. Finally, Chapter V summarizes the results and 
conclusions focusing on SAR applicability to ship detection. 
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11. BACKGROUND 


A. THE SYNTHETIC APERTURE RADAR (SAR) FUNDAMENTALS 

1. Basics 

SAR systems. 

A Synthetic Aperture Radar (SAR) is a Pulse-Doppler imaging radar that 
transmits short radio pulses in the microwave region of the spectrum. A SAR obtains 
resolution in range through travel time measurement, and uses Doppler shift of the 
backscattered signal to obtain resolution in azimuth as shown in Figure 1. The azimuth 
resolution is obtained using the concept of the synthetic antenna or synthetic aperture 
(Figure 2). A very long antenna is synthesized by moving a small antenna along a 
convenient path on a platform flight path or azimuth direction. The returning Doppler- 
shifted backscattered signal (amplitude and phase) is digitally processed to obtain 
information in the azimuth direction. The combination of time-resolved range and 
frequency-resolved azimuth results in a SAR image of target backscatter. 

The SAR systems create a large amount of data necessitating extensive processing 
to produce images with the required resolution. This process is usually performed and 
stored in ground stations. In the early days of SAR observations, the processing was done 
optically, but due to the huge increase in computational power, all processing is currently 
done digitally. 

The SAR technique was extended to use two antennas leading to the generation of 
three-dimensional images of the illuminated surface (Graham, 1974). This technique is 
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usually referred to as across track Interferometric Synthetic Aperture Radar (IFAR); 


however, IFAR is not the subject of this thesis. 



Figure I. Azimuth and Range direction. 

[From http://southport.jpl.nasa.gov/desc/imagingradarv3.html1 


' Syilhelc Aperture Length, L ~ 



Figure 2. Azimuthal Synthetic Aperture 

[From http://southport.ipl.nasa.gov/desc/imagingradarv3.html 
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2. SAR theory 

The foundation of SAR theory is the synthesized antenna aperture to increase the 
along track or azimuthal resolution. The SAR theory is described in Curlander and 
McDonough (1992), and Elachi (1982). Figure 3 in three dimensions and Figure 4 in two 
dimensions summarize all terminology and geometry for SAR. A full, detailed 
description of SAR theory is beyond the scope of this work but some basics are needed to 
better understand the ship detection problem described below. 
a) Azimuthal resolution 

The well-known result shown below is described by Curlander and 
McDonough (1992) among others and began with the azimuth resolution for Side- 
Fooking Aperture Radar (SEAR). The resolution is independent of the range and the 
radar wavelength, and improves with a smaller aperture, given by 



where da is the antenna real aperture in azimuthal direction. 

The corresponding resolution expression for an orbital SAR returns the complete 

result 


Pa 


Re + h 2 


[ 2 ] 


where Re is the earth radius and h is the platform altitude (Oliver, 1998). 
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b) Range resolution 

The minimum separation between two resolvable points in SAR can be 
defined as the minimum range separation of two points that can be distinguished by the 
system as separate. For SAR sensors 


n - 


cx 


sinr] 2 sinr] 


[3] 


or after pulse compression techniques 




26^ sin T] 


[4] 


where c is the light speed, is the resolution of the slant range, r] is the incident angle 
between the radar beam and normal to the earth’s surface and x ^ is pulse duration. 

c) Design considerations (constraints and limitations) 

Below is a summary of the highlights of the critical design aspects of a 
SAR as an emitting and receiving sensor. Other technical construction constraints related 
to SARs are beyond the scope of this work. 

(1) Choice of Pulse Repetition Frequency (PRF) or Pulse 

constraints. 


The correct pulse duration is needed to achieve a reasonable range 

resolution and a sufficient echo signal to noise ratio (SNR). Chirp pulses and 

compression techniques are used to achieve the optimum PRF. This is a severe constraint 

in designing practical SARs due to the compromise between improving azimuth 

resolutbn and providing a wider swath width. Therefore, ignoring the earth’s effect and 

orbit curvature for simplicity, we have from Robinson, (1985) the following: 
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c 


[5] 


IV y 

— =-< PRF < - 

2.VF.sin0 


where PRF is Pulse Repetition Frequency, W is the swath width, V is the platform 
velocity, ? is the angle between the slant direction and vertical (incident angle). The 
minimum PRF limitation is easier to visualize in terms of antenna theory because in those 
terms what determines the minimum PRF is the distance between successive array 
elements in the real antenna. If this distance is greater the X/2, the so-called grating lobes 
(replications of the arrays’ main lobe at increasing intervals on either side of the main 
lobe) will be produced. 

(2) Sidelobes Minimization. 

Stimson, 1998 states. 


Performance of a synthetic array radar may be degraded by both range sidelobes 
due to pulse compression and the sidelobes of the synthetic array. The sidelobes 
affect the radar maps in two different ways. First, the peaks of the stronger 
sidelobes may cause a string of progressively weaker false targets to appear on 
either side of a strong target. Second, the combined power of all sidelobes-called 
integrated sidelobe return-together with noise tends to fog or wash out the detail 
of the maps. 


(3) Motion Compensation 

Platform velocity and orbit are not perfectly constant as SAR 
theory assumes. The whole SAR concept is based on very slight signal phase changes 
over comparatively long periods. Therefore, it is essential that any motion during 
short periods be accounted for. A complete study of this error source is beyond the 
goal of this thesis, but the reader is referred to AGARD-LS-182 (1992). 
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(4) 


Antenna size limitations 


The reduction in antenna size is limited by the operational 
necessity to generate a sufficiently powerful signal in order to obtain enough returned 
signal. A minimum value of the vertical antenna dimension is needed to focus the 
beam into the desired swath width. These two design constraints are other practical 
limitations to SAR system resolution improvements. 

3. SAR operating modes 

There are three operational modes of SAR systems: STRIPMAP, SCAN AND 
SPOTLIGHT, pictorially sketched in Figure 5 (Franceschetti, 1999). These STRIPMAP 
is the most popular which was developed and adapted into the following modes: 

• Squint forward and backward - an array that, although it loses some 
azimuthal resolution, can execute applications like targeting/detection 
ahead of and behind the platform, and backscattering analysis of the 
properties of the illuminated surfaces with respect to the azimuth angle. 

• Doppler beam sharpening provides high quality, continuously updating 
maps of large expansions of the surface with the penalty of reducing cross¬ 
range resolution distance Multilook mapping. This mode improves map 
quality by averaging out the scintillation of the radar returns. 

• The SCAN mode tremendously increases the range swath dimension at the 
expense of azimuthal resolution while stepping the antenna beam to 
neighboring sub-swaths in range direction. 
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Figure 3. 
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SAR terminology and geometry 1 [From AGARD-LS-182,1992] 
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Figure 4. SAR terminology and geometry 2 [From AGARD-LS-182, 1992] 
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• The SPOTLIGHT mode is steered during overall acquisition time to 
illuminate the same area. A huge gain in azimuthal resolution is traded for 
a loss in coverage (small illuminated area along the sensor path). 





omtcvo*' 
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Figure 5. SAR operational modes [From Franchescetti, 1999] 
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B. 


SAR IMAGE PROCESSING 


Using Stewart,(1985), the SAR image proeessing produees high spatial resolutbn 
in the azimuthal direction and high spatial resolution in range. Processing the radar signal 
in essence consists of taking a scattering time history at each range and convolving it 
with the point-target response to produce a line in the image. Using a continuous 
approach, Franceschetti (1999) states the two-dimensional processed SAR image 
expression as follows: 



■') = ||y(x.r)sin: 


^ {x-x) 
Ax 


sin c 


f(r’-r) 

Ar 


dxdr 


[ 6 ] 


A 

where y(x, r) represents a two dimensional reflectivity scene pattern including the phase 
factor. Symbols x and r are the azimuthal and range direction respectively. 

I. Geometric Resolutions 

From Franceschetti (1999), also. 


Simply speaking, geometric resolution is the ability of the system to locali z e 
nearby objects. More precisely, the resolution length is the minimum spacing 
between two objects that are detected as separate entities, and are therefore 
resolved. 


Both results for range and azimuthal were already presented in [1] and [3]. 
a) Focused x Unfocused 

An “unfocused” SAR processing system is one which attains its along- 
track resolution by simple frequency filtering of the Doppler waveform. Higher azimuthal 
resolutions are only possible in a focused SAR system, which can account for the 
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variation in Doppler frequency of a target as it passes through the footprint resulting in 
the along-track resolution limit of half the antenna length. 

2. Radiometric Resolutions 

Franceschetti, (1999) defines radiometric resolution as “ measure of the ability of 
the system to discriminate, or resolve, areas of different scattering properties. These are 
described by the reflectivity pattern (y (x, r)) of the illuminated surface.” 


Franceschetti, (1999) also states that changes in the reflective pattern are related 
to two different processes, one macroscopic (image intensity modulation proportional to 
y{x,rY due to surface shape and/or its electromagnetic parameters as a function of space 
coordinates) and the other microscopic (essentially phase change due to surface 
roughness). A conventional definition of radiometric resolution given by Brooks and 
Miller, (1979) and Franceschetti, (1999) covers both microscopic and macroscopic 
process follows: 


^rd =10* log 


V 


[7] 


where p, is the mean value and a is the standard deviation of the image intensity 
referring to a distributed target of constant reflectivity. This states that adjoined areas 
with different intensities can be resolved, provided that their difference is larger than 
Ard. 


3. Radiometric and Geometric Distortions 

In practice, many factors can affect image quality by causing signal amplitude or 
phase modulation. The former is a consequence of antenna failure in stabilization 
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creating a change in the antenna point angle. Phase modulation is a consequence of an 
unpredicted change in the range between the object and the antenna, giving rise to an 
unknown receiving signal, which affects azimuthal pulse compression. 

Oliver, (1998) summarizes that the amplitude modulation leads to radiometric 
distortions and increased side lobes, and phase modulation leads to image defocus, 
geometrical distortions (foreshortening, layover and shadow), radiometric distortions and 
increased sidelobe levels. Image defocus, radiometric distortion, and increased sidelobe 
levels restrict target detection and recognition capability. Geometrical distortions prevent 
direct pixel-to-pixel comparison, which is relevant to multichannel registration, target 
change detection, and direct image comparison. Motion Compensation (MOCO), made 
by accurate accelerometers and inertial navigational units mounted on the antenna, can 
reduce amplitude modulation effects. However, after this processes some residual 
unknown phase modulation still exists. 

4. SAR Signal Statistics 

In the case of a SAR scene, the only way to describe its roughness and density of 
scatters, for surface and volumetric scattering respectively, is in terms of statistical 
parameters. Data characterization throughout a probability density function (PDF) is a 
matter of compromise between simple models with few degrees of freedom and complex 
models with many degrees of freedom. Some well-known PDF with two-degrees of 
freedom like Weibull, lognormal, and the K distribution are used to statistically represent 
a SAR image. We will focus on the K-distribution since it is the PDF used in our ship 
detection model. 
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5. Speckle 

The fluctuation in the receiving signal, known as fading, caused by the diverse 
nature of the return echo, generates a noise like grainy appearance on SAR imagery 
called speckle. It influences the ability to estimate image properties and thus is central to 
information retrieval from individual SAR images. Speckle has a noise like appearance, 
but it is not noise. An electromagnetic measurement can be useful, for example, in SAR 
interferometry. It can also be understood as an interference phenomenon in which the 
principal source of the noise like quality of the observed data is the distribution of the 
phase terms. 

C. RELEVANT TOPICS IN SHIP DETECTION 

I. Radar Cross Section (RCS) 

The RCS describes the backscattering property of the target and depends on its 
size, shape and orientation as well as on wavelength and polarization of the incident 
signal (Skolnik, 2001). From the definition of RCS (a) the cross section of an object that 
scatters equally in all directions (isotropic) is equal to its projected area. 

a) Normalized Radar Cross Section (NRCS) 

The quantity o°. Radar Cross Section per Unit Area, provides a 
normalized parameter that can be used to describe RCS of the Sea independent of the 
radar illumination. 

The normalized radar cross section (NRCS) is given by 

o^^o/A [7] 

where A is the area of a smooth surface that corresponds to the mean land or sea surface 


area contained within the radar’s cell resolution. The illuminated area is a complicated 
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function of radar system parameters and is generally predictable as is described below in 
Section 3. 

2. Ship Radar Cross Section 

Ship reflectance is a complicated problem without a general solution. Skolnik 
(2001) proposed an empirical relation between displacement in tons and ship radar cross 
section. Vachon et al. (1997) using the Skolnik relationship developed the following 
linear regression between ship weight and length for the Bedford Institute of 
Oceanography Fleet and for some ships participating in the MARCOT’95: 

a = D = O.OS.L^'^ [8] 

where L is ship length and D is displacement. 

Olsen et al. (2000) improved the above formula including the effect of incident 
angle using the same data from MARCOT’95 

a = D = R(0).O.O8.L^^^ [9] 

where 


R(0) = 0.78+0.11.0 [10] 

All these formulations show that ship RCS is founded unobserved values and empirical 
methods. 

3. Ocean Radar Cross Section 

a) Nature ofo° for Ocean and Incident Angle Dependency 

The NRCS is sensitive to the character of the sea surface and is highly 

variable. It can change as much as 10 dB in a one-minute interval due to variations 


among measuring instruments and their calibration Long (2001). Therefore, measurement 
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errors and uncertainties tend to obscure the weak functional relationship between the 
wavelength (^) and o° for most incident angles. This relationship is a complicated 
function of incidence angle. For instance, at low incidence angle near the vertical, o° is 
highly dependent on sea state and, for high incident angles, o° is a strong function of X . 
For horizontal polarization, the curve of o° versus the grazing angle ^ 0 for low 
frequencies is expected to be the same for microwaves except for a larger critical angle. 
Therefore, the curve for horizontal polarization might appear as in Figure 6. 

b) Dependence of <3° on Polarization 

Polarization affects a° for different surface conditions or for different 
wavelengths. From Long (2001) makes the following statements that hold for the ratio 
between a vv (NRCS for vertical transmitted and vertical received polarization) and a hh 
(NRCS for horizontal transmitted and horizontal received polarization): 

1. The ratio increases with an increase in wavelength, 

2. The ratio decreases with an increase in sea roughness, and 

3. o°hh can exceed o'w for heavy seas and small depression angles. 

It should not be erroneously assumed from the above that under very rough 
conditions o w and o hh are about equal for all incident angles. 

c) Dependence of o° on Wind and Sea 

Some of the sea surface characteristics known to influence various 
features of radar echos include the period and shape of the waves, the wave height, wind 
ripples, and the presence or absence of whitecaps and spray. However, observations 


^ 0° incidence angle from nadir (vertical) is equal to 90° grazing angle measured from horizontal plane 
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indicate that average and medium values of o° inerease as the sea beeomes rougher for 
all grazing angles except for those near vertical incidence. Vaehon et al. (1997) uses an 
empirical polarization ratio constructed by Unal et al. (1991) to transform from modeled 
values of o°vv to values of o°hh neeessary for the ship length eomputation. 



Figure 6. Shape of a° versus Grazing Angle curves for horizontal polarization 
[From Long, 2001] 
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m. SHIP DETECTION PROBLEM 


A. SHIP DETECTION APPROACHES 

I. Ship-Generate Surface Waves (Wake detection) 

Ship wake patterns exhibit various characteristics in SAR images under diverse 
environmental conditions. Wake structures in an image can be classified in three 
categories. First, surface waves generated by the ship; second, turbulent wakes or vortices 
containing relatively persistent but non-propagating currents visible in SAR imagery 
through the interaction of ambient waves with these surface currents (Lynden, 1988). 
Third, internal waves generated by the ship that are visible due to interaction with the 
short surface waves. 

Ship wakes that are observed in SAR imagery are created by two mechanisms. 
First, variations of short waves (centimeters-scale) appearing as bright narrow “V” wakes 
are directly observed due to Bragg scattering in images; second, longer waves (decameter 
scale) forming the classical Kelvin-wake system are observed as they modulate the 
formation of Bragg scattering waves. Standard ship wake detection techniques exploit 
the first subcategory due to better pixel intensity (bright) response in SAR images. Some 
techniques (Rey et al. (1990), Copeland et al. (1995) and others) have been developed to 
use the Radon transform^ in extracting linear features that are correlated with ship wakes. 

Many specialists consider point target approaches a more efficient “pattern 
recognition” technique than wake detection. Wake detection techniques are very useful, 
after a candidate vessel is chosen, in order to confirm or reject the target as ship. 

^ A commonly used technique for finding lines in an image. Very robust in the presence of high levels 
of noise but lacking in differencing between short and long lines or end points of short lines. 
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Furthermore, in key applications like fisheries monitoring the vessels may often travel 
too slowly to generate significant wake signature (Olsen, 2000). This work will present a 
performance analysis of new spacebome SARs using a point target detection approach. 

2. Point Target Detection (Ship detection) 

The main task of point target detection is to detect a single or a small cluster of 
pixels with the higher backscatter intensity from a ship within the background scene 
intensity. Generally, the background ocean has a lower backscatter intensity than ships. 
The ocean backscatter varies significantly with wind speed and wind direction relative to 
the SAR view angle, SAR beam incidence angle, and polarization of the returned power. 
The challenge of ship detection is unambiguously identifying the small number of bright 
ship pixels among the variable intensity of the background population of pixels. 

3. Ship Detection Techniques 

An intuitive approach to ship detection is to divide the scene into small frames 
assuming a single ocean backscatter value within each frame that determines a detection 
threshold for candidate vessels. The simplest method is to set the threshold intensity at n 
standard deviations above the image mean intensity. One variation of this technique uses 
two windows inside the image with different sizes. The windows are moved across the 
entire image using the larger window to compute the ocean background statistics while 
the smaller window is used for the actual target search. 

Some more sophisticated parametric approaches fit a PDF (usually as a K- 
distribution ) to the data and estimate the threshold based on the uncertainty level desired. 
One non-parametric method uses a probabilistic Neural Network approach. All these 
methods are designed to reach a balance between detection rate and the occurrence of 
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false alarms. To achieve this, some analysis methods are applied such as contrast 
measures, homogeneity tests, morphological filters, wake analyses, and combinations of 
these. 

a) The Constant False Alarm (CFAR) method 

The CFAR method sets a threshold for radars operating in a sea-clutter 
environment based on a PDF for the image ocean background intensities, integrated until 
a required significant level is reached. This significant level, Tjc, corresponds to a CFAR 
of l-Tjc for the intensity PDF 

pix)dx [8] 

Two main problems arise from this method. The first problem is to find a PDF that fits an 
image intensity histogram. Second, F must be computed. Vachon (1997) uses Tjc = 
0.995 and p(x) is defined as a K-distribution for various order parameters. 
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IV. NEXT GENERATION SHIP DETECTION EXPECTED 

RESULTS 

The main question to be answered by a Ship Detection Model is what size of 
vessels can be detected in SAR imagery and how reliable is the detection under various 
observational conditions (Olsen, 2000). Secondarily, the model can indicate whether 
information like type, speed and heading can be retrieved from SAR imagery. 

The key operational requirements for the next generation spacebome SARs that 
are needed to improve ship detection are: 

• Increase in resolution to be applied in target detection; 

• increase in swath coverage to allow widespread ship monitoring and surveillance; 

• larger incident angles (> 35°) to reduce ocean clutter and increase ship 
backscattering; 

• greater coverage frequency in combination with higher resolution and 
programmable modes for selected areas. 

Almost all of these desirable aspects will be flying on the next spacebome SARs. The 
goal of this thesis is to apply, in a reliable model, the new specifications, mainly for 
resolution and incident angle, and retrieve the general performance behavior for different 
satellites and their different modes. 
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A. 


RELEVANT NEXT GENERATION SATELLITE PARAMETERS 


This section gives a brief description of all three satellite and sensor 
characteristics and the launch schedule over the next 2-3 year window. Details for 
RADARSAT2, ENVISAT/ERS2 and ALOS/PALSAR SAR systems are summarized in 
the following sections 

1. RADARSAT- 2 

Table 2 presents RADARSAT-2 values for swath width, incident angles, Number of 
Looks, azimuthal and range spatial resolutions, multipolarization capability and others 
characteristics, like orbit parameters and coverage access using 500km swath width. 
Figure 7 shows the coverage pattern of all operational modes and Table 3 shows 
RADARSAT-2 inovations. 



Figure 7. RADARSAR 2 Modes [From CCRS/RADARSAT- 2 site] 
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Table 2. RADARSAR 2 characteristics [From CCRS/RADARSAT- 2 site] 


Standard Beam 
Modes 

Nominal 

Swath 

Width 

Incidence 

Angles 

Number 
of Looks 

Approx. 

Resolution 

Wide 

100km 

20-50 

1x4 

25m X 28m 

Eow Incidence 

150km 

20-45 

1x4 

25m X 28m 

High Incidence 

170km 

10-20 

1x4 

40m X 28m 

Eine 

70km 

50-60 

1x4 

20m X 28m 

ScanSAR Wide 

50km 

37-48 

1x1 

10m X 9m 

ScanSAR Narrow 

500km 

20-50 

4x2 

100m X 

100m 

Standard Quad Pol. 

300km 

20-46 

2x2 

50m X 50m 

Eine Quad Pol. 

25km 

20-41 

1x4 

25m X 28m 

Triple Pine 

25km 

30-41 

1 

1 Im X 9m 

Ultra- fine Wide 

50km 

30-50 

3x1 

1 Im X 9m 

Ultra-fine Narrow 

20km 

30-40 

1 

3m X 3m 


ORBIT CHARACTERISTICS 


Altitude (average) 

798 km 

Inclination 

98.6° 

Period 

100.7 minutes 

Ascending Node 

18:00 hrs 

Sun-synchronous 

14 orbits per day 

Repeat Cycle 

24 days 


COVERAGE ACCESS USING 500km SWATH WIDTH 


North of 70° H 

Daily 

North of 48° H 

Every 1-2 days 

Equator 

Every 2-3 days 
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Table 3. RADARSAT-2 innovations [From Sanden (2001)] 


INNOVATIONS 

BENEFITS 

GPS receivers onboard 

+60-meter real-time position 

information 

10 ms delay between imaging 
modes 

Faster mode changes 

Yaw-steering for zero-Doppler 
shift at beam center 

Facilitates image processing 

Higher downlink power density 

3-metre minimum size antenna on 
ground allowing station portability 
and Lower "cost of entry" for new 
ground stations 

3-meter ultra-fine resolution 

Highest-resolution commercially 

available SAR 

Left-and right-looking capability 

Faster revisit time, 2000 km 
accessibility swath and Routine 
Antarctic mapping available 

Fully polarimetric imaging modes 

Enhanced capabilities for various 
applications 

Solid-state recorders for onboard 
image storage 

Higher reliability, faster image access 
and Simultaneous record and 
downlink 


2. ENVISAT/ASAR 

ASAR has five mutually exclusive modes of operation which could be classified 
in two categories: Global mission subdivided in to Global monitoring and Wave mode 
that have a low data rate and hence an operation capability up to 100% of the orbit, and 
regional mission subdivided in narrow swath modes, image mode and alternating 
polarization mode. Table 4 summarizes some operational parameters and Table 5 the 
mode’s characteristics. 
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Table 4. 


ENVISAT/ASAR Operational Parameters [From ENVISAT Internet site] 


Instrument 

Parameters 

Image Mode 

Aiternating 

Poiarization 

Wide Swat 
h 

Giobai 

Monitoring 

Wave Mode 

Swath width 

Up to 100 km 

Up to 100 
km 

> 400 km 

> 400 km 

5 km vignette 

Operation 

time 

Up to 30 min per orbit 

Rest of orbit 

Data Rate 

Up to 100 Mbit/s 

0.9 Mbit/s 

Power 

1365 W 

1395 W 

1200 W 

713 W 

647 W 


Table 5. ENVISAT/ASAR Modes [From ENVISAT Internet site] 


ASAR 

Swathes 

Range 

Incidence 

Angle 

Swath Far 

Width Range 
[km] 

Near Incidence 
Angle 

IS1 

108.4 - 109.0 

14.1 - 14.4 

22.2 - 22.3 

IS2 

107.1 - 107.7 

18.4 - 18.7 

26.1 -26.2 

IS3 

83.9 - 84.3 

25.6 -25.9 

31.1 -31.3 

IS4 

90.1 -90.6 

30.6 -30.9 

36.1 -36.2 

IS5 

65.7 -66.0 

35.5 -35.8 

39.2-39.4 

IS6 

72.3 - 72.7 

38.8 -39.1 

42.6-42.8 

IS7 

57.8 -58.0 

42.2 - 42.6 

45.1 -45.3 


3. ALOS/PALSAR 

The Advanced Eand Observing Satellite with a Phased Array type E-band 
Synthetic Aperture Radar (PAESAR) is the Japanese next generation spacebome SAR. 
PAESAR uses Eband frequencies, different from RADARSAT and ENVISAT/ASAR, 
and has land observation as a major task and therefore has a variable beam elevation. The 
development of the PAESAR is a joint project between NASDA and Japan Resources 
Observation System Organization (JAROS). Its major specifications are listed in Table 6. 
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Table 6. ALOS/PALSAR specifications 


[From ALOS/PALSAR site] 


Observation Mode 

Eine Resolution Mode 

ScanSAR Mode 

Frequency 

E-band 

Polarization 

HH or VV (option: HV or VH) 

Spatial Resolution 

10m (21ooks)/20m (4 looks) 

100m 

Swath Width 

70km 

250-360km(3 - Sscans) 

Off- nadir Angle 

18 - 48deg. 

S/A 

25dB 


-25dB 


B. SHIP DETECTION MODEL APPLIED ON RADARSAT-I 

Vachon et al. (1997) develops a model to study RADARSAT- 1 (C band HH 
polarization) performance for ship detection using the various beam modes. The model 
was validated with data acquired during the RADARSAT-1 ship detection/validation 
field program. Minimum detectable ship length is used as a Figure of Merit (FOM) for 
ship detection using RADARSAT-1 standard modes and ERS-1 modes. Validation of the 
model’s key assumptions are provided, including those concerned with the hybrid C-band 
HH polarization ocean radar cross section model, the image probability density function 
and the ship radar signature model. This model is being updated by the Canadian Center 
of Remote Sensing (CCRS) and was applied and tested in the Ocean Monitoring 
Workstation - a system developed by Satlantic, Inc. with technical and financial support 
of CCRS, The Department of Fisheries and Oceans (DFO), Canadian Coast Guard 
(CCG), the Department of National Defense (DND) and the Canadian Space Agency 
(CSA). The OMW is designed to provide operational users of marine data with near-real 
time, value added ocean information derived from RADARSAT-1 SAR images (Vachon, 
2000 ). 
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1. Model Assumptions 

The model is based on three assumptions: 

(1) the SAR image probability density function for ocean scenes can be 
described by the compound K-distribution model for sea clutter; 

(2) the radar cross section of the ocean can be estimated given wind speed and 
beam geometry of C-Band radars; 

(3) ship size is directly related to the radar cross section of the ship. 

The model was developed using ERS-1 SAR data prior to RADARSAT-1 launch and 
was subsequently validated with RADARSAT-1 data. 

2. SAR Ocean Scene Image Statistics 

The compound K-distribution model is an empirical model for sea clutter that is 
apphcable to high-resolution radars that resolve fine structure on the sea surface. 
Gaussian distributions do not fit for this kind of processes and the K-distribution does so 
due to the following observations from Rey, 1996: 

The envelope time history of a demodulated signal from an individual range cell 
shows fast fluctuations from pulse to pulse that is modulated by an underlying 
structure. 

The fast fluctuation component, or ‘speckle’, decorrelates from pulse to pulse 
with frequency agile signals but is correlated for between 5 and 10 ms with fixed 
frequency signals. This implies a large number of scatters within each illuminated 
patch. By applying the central limit theorem, the ‘speckle’ can be approximated 
by the Rayleigh distribution. 

The Slower modulation component is the local sea clutter mean level, which has a 
long temporal decorrelation period and is not affected by frequency agility. It 
characterizes the mean level variation of clutter, including clutter ‘spikes’, and its 
intensity is found to be gamma distributed; 

The gamma distribution is found in this model to be the first term in a Laguerre 
series expansion of the modulation PDF. The implication is that the non-Gaussian 
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nature of sea clutter seems to be due to more to scatter bunching by the sea wave 
structure than anything else. 


The proposed distribution to describe the SAR PDF of image intensity is the 
multi-look intensity K-distribution (Oliver, 1993): 


p{l)=] 



(z,+d)/2 

V 

9 

( LvI 3 

1/2 

(</>J 

r(L)r('u) 


[</>J 



[9] 


where p(I) is the probability of the image intensity 1, <1> is the mean image intensity, L 
is the number of statistically independent looks, V is an order parameter for the intensity 
modulation process, F is the gamma-function, Kv-l is the modified Bessel-function of V-L 
order. 


Using the above PDF, we can now compute the critical intensity using [8] for any 
significant level, Tjc- To retrieve the same results as Vachon et al. (1991) we use 0.995 
significance level, a trapezoidal numerical routine for integration (0.005 and 0.01 
intervals), and image mean </> = 1. 

3. RADARSAT-1 SAR parameters 

The RADARSAT-1 parameters are summarized on Table 7. The values are for 
nominal (near beam center) parameters for the 21 modes of the RADARSAT-1 SAR. 
ERS-1 values are also given for comparison. 

4. Radar Cross Section of the Ocean 

A physically-based model for the a° relevant to RADARSAT-1 C-Band HH 
polarization is not available. Vachon et al. (1997), using Gower et al. (1993), based on 
Snoeij et al. (1991), incorporates empirical results from Unal et al. (1991) for backscatter 


32 



polarization ratio (PR“). Table 8 lists the conversion factors to convert from a VV to HH 
Ocean Radar Cross Section polarization model. 

Table 7. RADARSAR-1 parameters [From Vachon (1997)] 


Beam Mode 

Swath [km] 

Looks(L) 

0[deg] 

Padm] 

PrJm] 

o°„ddB] 

Dchirp [km] 

SI 

100 

3.1 

23.5 

27 

24.2 

-23 

15.8 

S2 

100 

3.1 

27.5 

27 

20.4 

-21 

13.6 

S3 

100 

3.1 

33.5 

27 

25.3 

-23 

11.4 

S4 

100 

3.1 

37.0 

27 

23.4 

-23 

10.5 

S5 

100 

3.1 

39.0 

27 

22.1 

-22 

10.0 

S6 

100 

3.1 

43.5 

27 

20.3 

-24 

9.1 

S7 

100 

3.1 

47.0 

27 

19.1 

-23 

8.6 

FI 

50 

1 

38.5 

8.4 

8.3 

-23 

10.1 

E2 

50 

1 

40.5 

8.4 

7.9 

-22 

9.7 

E3 

50 

1 

42.5 

8.4 

7.6 

-23 

9.3 

F4 

50 

1 

44.5 

8.4 

7.3 

-24 

9.0 

F5 

50 

1 

46.0 

8.4 

7.1 

-23 

8.6 

W1 

150 

3.1 

25.5 

27 

33.8 

-22 

14.6 

W2 

150 

3.1 

35.0 

27 

24.6 

-22 

11.0 

W3 

150 

3.1 

42.0 

27 

20.8 

-25 

9.4 

EHI 

75 

3.1 

50.5 

27 

18 

-25 

8.1 

FH2 

75 

3.1 

51.5 

27 

17.7 

-25 

8.0 

EH3 

75 

3.1 

53.5 

27 

17.3 

-25 

7.8 

EH4 

75 

3.1 

55.5 

27 

16.8 

-23 

7.6 

FH5 

75 

3.1 

57.0 

27 

16.6 

-23 

7.5 

FH6 

75 

3.1 

58.5 

27 

16.4 

-23 

7.3 

ELI 

170 

3.1 

16.5 

27 

39.1 

-22 

22.2 

ERS-I 

100 

4.9 

23.0 

30.7 

38.4 

-24 

14.2 


This table includes the following parameters and modes: Sstandard; Ffine; W- wide; EH- 
extended high incidence; EL-extended low incidence; L is the number of independent looks; 0 is the 
incidence angle paz is the azimuth resolution; prg is the ground-range resolution; a°ne is the noise- 
equivalent clutter level; Dchirp is the ground-projected chirp length. 


Table 8. Vertical to Horizontal C-Band (5.3 GHz) Polarization Ratios [From Unal 

etal. (1991)] 


Inc.ang/ 

Wind(m/s) 

20° 

30° 

45 

2 

1.05 

2.07 

4.78 

4 

0.79 

2.37 

5.25 

6 

0.65 

2.57 

5.47 

8 

0.56 

2.70 

5.59 

10 

0.51 

2.80 

5.66 

12 

0.49 

2.88 

5.71 

14 

0.49 

2.95 

5.75 
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The well-developed C-band VV (vertical/vertical) polarization model (o°vv) for 
ERS-1 scatterometer wind retrieval (the C band Models-CMOD series), was used here in 
conjunction with Table 8 ratios to produce the HH (horizontal/horizontal) RCS of the 
Ocean (o°hh) as showed below 

o°hh = o°vv-PR“ [10] 

5. Radar Cross Section of Ships 

Based on Skolnik (2001), Vachon et al. did a linear regression between ship 
weight and length for the Bedford Institute of Oceanography (BIO) fleet and some ships 
participating in the military exercise MARCOT’95, resulting in the following empirical 
relation: 

aship= = O.OS.L^^^ [11] 

where Oship is the ship radar cross section, ‘D’ is the ship displacement and L’ ship 
length. 

The model can then determine the ship detectability by a SAR based on a ship 
scale parameter. Knowing the critical intensity (T) of the relevant PDF from [9] and the 
ocean normalized cross section (a°) from [10], the model can determine the 
corresponding minimum point target RCS for detection at a chosen probability level as 

[ 12 ] 

where Paz is the satellite azimuthal resolution cell size and Prg is the ground range 
resolution cell size. 
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6. Ship Length computation 

Combining [11] and [12] we have the minimum detectable ship length as a 
function of the Ocean RCS, resolutions and critical intensity, as follows: 

L = (I^c>V.,P.,/0.08)''' [13] 

Applying this results to RADARSAT-1 data, Vachon et al. (1997) produced results for 
SI, S3 and S7 RADARSAT-1 modes and ERS to test the model, considering -23 dB as 
the RADARSAT-1 noise equivalent normali z ed radar cross section ^ne) with clutter 
noise ratio (CNR) = o°/One and One = -24 dB for ERS. 

7. Vachon et al. Ship Detection Model Results 

Some results for this model for RADARSAR 1 mode SI, S3, S7 and ERS SAR 
are presented in Eigure 8. Each plot shows the minimum detectable ship length as a 
function of wind speed using equation [13] where solid lines represent the Upwind 
direction (ie, (|)=0, the largest a°, hence the worst case - least detectable due to a bright 
background) while dashed lines represent the Crosswind direction (ie,(|)=90°, the smallest 
o°, hence the best case). The three curves represent variation of the order parameters for 
K-distribution from top to bottom: v= 4, 10, 8 (ie, from worst to best case). To facilitate 
comparison among the various RADARSAT-1 beams and modes, Vachon et al. (1997) 
defines a EOM for ship detection using Upwind (|) = 0° direction, wind velocity U = 10 
m/s, and order parameter v = 4. The resulting plot is shown in Eigure 9. 
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8. 


Vachon et al. results validation 


A validation field program was held in March/April 1996 to provide quantitative 
information to validate the ship detection model. Some of the elements used to validate 
the model are showed in Vachon et al. (1997) and summarized here. They include 
acquisition of wind and wave data from two buoys, acquisition of nine RADARSAT-1 
SAR scenes in various modes, placement of ship near the buoys at the time of 
RADARSAT-1 SAR passes, and surveillance flights over the area to identify additional 
ship targets on an opportunity basis. 


After RADARSAT-1 calibration and SAR data processing, Vachon et al. (1997) 
concludes. 


We validated some of the ship detection model assumptions using our calibrated 
RADARSAR SAR data. In particular, we show that there is excellent agreement 
between the observed RADARSAT-1 SAR cross section and those predicated by 
the hybrid model. Furthermore, the K-distribution in shown to be a suitable PDF 
for RADARSAT-1 SAR ocean images. We also show that our very simple ship 
length-dependent cross section model is within the correct order of magnitude 

C. SHIP DETECTION MODEL APPLIED TO NEXT GENERATION 

SPACEBORNE SAR PARAMETERS 


I. Some Considerations 

From Vachon et al. (1997), we have some important statements about the radar 
cross section of the ship. They are: “...For the ship considered, we find that our simple 
relationship is within 20dB, but that the anomaly tends to increase with increasing 
incidence angle”, and more “ ...Although we have not considered ship orientation at all 
in this analysis, use of the best fit line to our data allows recovery of the vessel length to 
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Figure 9. Figure of Merit from Vachon et al. (1997) 

within 20% based on the measured ship cross section alone. It is apparent that our simple 
ship model will tend to underestimate the ship detectability.” 

Due to the good results presented and supported by Vachon et al. (1997), the 
model is used here to describe the physical processes that limit ship detection and to 
quantify the ship detection improvements that can be expected in the next generation 
spacebome SARs. 

Unal’s et al. (1991) table is used as reference and including only values linearly 
interpolated inside the table wind range (except for mode S7 where an linear 
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extrapolation to 47° of incident angle) The computations begin at 2 m/s wind velocity for 
to two reasons. First, due to a lack of Unal et al. (1991), values inside this interval to 
confirm our values. Secondly, because our computed curves have a better agreement with 
Vachon et al. (1997), only results between the 2 to 14 wind speed interval are included 
here. 


2. Reproducing Results with Vachon et al. model: 

The Vachon et al. (1997) model is reproduced here using RADARSAT-2 
specifications and the same image parameters used for RADARSAT-1 except for 
resolution. All considerations taken in IV-C-1 were generated during this computation 
due to a small difference in length values and a more linear result for all RADARSAT-2 
and ERS modes. Figures 10,12,14 and 16 from below with Figure 8. 

3. Expected Result for all RADARSAT 1 modes kept in 
RADARSAT-2 

All four cases resulting from [13] show the same ship length curve behavior as 
wind speed increases for the same fixed values used by Vachon et al. (Figure 8). The 
unchanged beam parameters were Number of Fooks (3.1), each mode incident angle, and 
azimuthal and range resolution. Figures 11,13,15 and 17 show the cumulative K- 
distributions for each case that were used to produce an critical intensity value for each 
mode. 

Table 9 shows the comparison between Vachon et al. (1997) results and the 
results produced for this thesis. The differences in minimum detectable ship length are 
typically less than 5 m. Therefore the reproduction of the model is acceptable for the 
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main purpose of this thesis, the performance prediction for next generation ship detection 


spacebome SAR. 
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Figure 10. SI mode ‘Wind versus Ship’ length thesis results 
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Figure 11. SI mode Cumulative PDF with Critical Intensity Values 



Figure 12. S3 mode ‘Wind versus Ship’ length thesis results 
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Figure 13. S3 mode Cumulative PDF with Critical Intensity Values 





Figure 14. S7 mode ‘Wind versus Ship’ length thesis results 
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Figure 15. S7 mode Cumulative PDF with Critical Intensity Values 


Table 9. Vachon versus thesis length differences for SI, S3, S7 and ERS modes 
Obs.: Values rounded to meter and values computed for wind = 2 m/s 


Modes 

Sly 

Sit 

Slv-Slt 

S3y 

S3t 

S3v-S3t 

Worst 

27 

29 

+2 

13 

12 

-1 

Best 

20 

24 

-4 

10 

9 

-1 

Modes 

S7,. 

S7t 

S7y- S7t 

ERSv 

ERS, 

ERSv- ERS, 

Worst 

9 

4 

-5 

40 

38 

-2 

Best 

7 

3 

-4 

28 

30 

+2 


Worst is Upwind with Order Parameter=4; Best is Crosswind with Order Parameter = inf 
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Figure 16. ERS ‘Wind versus Ship’ length thesis results 
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Figure 17. ERS mode Cumulative PDF with Critical Intensity Values 
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4. Expected Results for Fine Mode Resolution 

Figures 18, 19 and 20 show the results using the RADARS AT-2 fine mode 
resolutions. The mean incident angle for each fine mode and the released Number of 
Looks, listed in Table 2, were used for each mode. All the Figures show the expected 
minimum ship length behavior as wind speed increases. Values with crosswind (best) 
and upwind (worst) conditions are shown with three different K-distribution order 
parameters. Despite the disadvantage of lack of coverage due to a small swath width, all 
expected detected ship length values associated with new high resolutions for fine modes 
(1.5 to 10 m) will lead not only to ship detection but also ship recognition. 
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Figure 18. Ship length for RADARSAR 2 Triple Fine mode 
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Figure 19. Ship length for RADARSAR 2 Ultra-Fine Wide mode 
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Figure 20. Ship length for RADARSAR 2 Ultra-Fine Narrow mode 
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5. Ship Length versiB Resolution Qualitative Expected Behavior 

Figure 21 shows how detectable ship length changes with a fixed wind velocity 
using one RADARSAR-2 mode (Standard), varying with resolution for 
upwind/crosswind (worst/best case), and with two Number of Looks values. Three major 
differences are noticeable. First, the high incidence angle detection can be increased in 
wide swath modes as much as five times compared to low incidence angle observations 
(ScanSAR narrow). Next, intermediate resolution modes (Standard), have less incident 
angle influence and almost the same influence of sea RCS and wind incident angle. 
Finally, at high resolution all three variables have about the same influence on final ship 
length computation. 

6. Ship Length versus Incident Angle Qualitative Expected Behavior 

Figure 22 shows the variation of detectable ship length versus incidence angle for 
a fixed wind velocity and RADARSAR mode (standard in this case) varying with 
resolution, upwind/crosswind (worst/best case) and Number of Looks. With an increasing 
incident angle the ship detectability increases by a factor of almost 2.5 for this specific 
case. High incident angle we have the best results for ship detection. 

Figure 21 and 22 confirm the major role of resolution and incident angle in the 
Vachon et al. (1997) ship detection model. The effect of resolution decreases as the 
incident angle increases. Number of Looks and sensor-relative wind direction have a 
secondary, but not disregardable influence. 
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Figure 21. Ship Length versus Resolution Results 
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Figure 22. Ship Length versus Incident Angle Results 
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7. Practical Application 

Given the wind condition as in Figure 23 at ON, 34W and 4N, 34W (North 
Brazilian Coast) from PIRATA project, we can use the Vachon et al. (1997) model to 
estimate ship detection performance. For instance, using the same parameters and 
statistics as in the tested model for the SI mode case on RADARS AT 2, we can expect 
that SAR images taken under the same conditions as in Figure 10 (mode wind velocity of 
6 m/s), ship lengths from 32 to 43 (best and worst case) are expected to be detected. This 
short example testifies to the model’s usefulness in analyzing real situations setting the 
best parameters for ship detection with RADARSAT 2. 




Figure 22. Wind histogram from PIRATA site 
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V. CONCLUSIONS 


The ship detection model outlined by Vachon et al. (1997) was used to verify and 
quantify ship detection improvements expected in the next generation spacebome SAR 
imagery. This thesis focuses on the influence of incidence angle, wind speed and 
resolution, on the detectability of ships in SAR imagery. The results are compared to 
RADARS AT 2 parameters to test the improvement expected from new capabilities. The 
impact of new polarized modes in ship detection is not tested in this thesis. 

Vachon et al. (1997) uses a model of the ocean Radar Cross Section (RCS) as a 
function of wind speed and incidence angle. The minimum ship length that can be 
detected against the background ocean radar cross section is determined by the critical 
intensity level obtained from a statistical relationship between ship size and RCS. 

The minimum detectable ship length is strongly dependent on wind speed due to 
the changes imposed by the wind on RCS of the ocean and on the critical intensity. For 
instance in RADARSAT-2 SI mode, with 3.1 Number of Looks, with 23.5° of incident 
angle with 99.5 % of confidence, we expect winds to change from 2 to 10 m/s while the 
detectable ship lengths change from 29m to 54m in the worst case and from 24m to 37m 
in the best case. For areas surveyed by the PIRATA project on the North Brazilian Coast 
(ON, 34W and 4N, 34W) with mode wind velocity of 6 m/s, minimum detectable ship 
lengths from 32 m to 43 m (best and worst cases) are expected to be detected. 

The minimum detectable ship length changes dramatically with beam incident 
angle. The increase in incidence angle reduces Ocean RCS due to a reduction of Bragg 
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scatter to the sensor and the increase of the ship RCS due to an increase in double bounce 
scatter. For instance, for RADARSAT 2 standard mode, with 10 m/s wind speed in the 
worst case (upwind related to the antenna), increasing the incident angle from 20° to 45° 
ehanges the deteetable ship length from 99m to 20m - nearly a faetor of five ehange. 

The minimum detectable ship length is strongly dependent on radar resolution 
since reduction of the pixel size increases the probability that photons from that pixel are 
reflected by the ship. For instance, changing modes in RADARS AT 2 from ScanSAR to 
UltrafineNarrow, thereby increasing resolution by a factor of 10 in area, with 10 m/s 
wind speed, 20° of incident angle and Number of Looks equal 1, the deteetable ship 
length will deerease from 280m to 15m (worst case) - almost a faetor 20 change. 

The best ship detection results are obtained with RADARSAT 2 Fine Resolution 
Modes. The minimum detectable ship lengths of 1.5 to 3.5 meters for fine mode 
resolution not only illustrate the improvement of ship detection, but also make ship 
recognition possible given the right conditions. Additionally, ship recognition will be 
improved by cross polarization. 

The general level of improvement for RADARSAT 2, ALOS/PALSAR and 
ENVI/ASAR arises from the new resolution modes and eross polarization Some specific 
design factors should be mandatory for future SARs designed specifically for ship 
detection. High incident angles orbits, shorter repeat cyele, multipolarization and large 
swath width at high resolution are among the main SAR characteristics needed for future 
ship deteetion satellites. 
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APPENDICES 


A. UNAL’S TABLE INTERPOLATION CODE 

function [Pr_Interpol_in_dB] = UnalsInt(V,theta,IntModeType) 

0 

% Interpolation to find Unal's et al. tabulated polarization ratio 

o 

% inputs 

% V - the wind speed in m/s 

% theta - the incidence angle in degrees 

% IntModeType - Matlab interpolation type, ie, 'linear','cubic',etc 

% output 

% Pr_Interpol_in_dB = Polarization Ration Interpolated in decibels 

Beam_Inc_Tb= [20 30 45 47]' 

Spd_Wind_Tb= [2 4 6 8 10 12 14] 

% Unal's et al., table values in dB extrapoleted until 47 degrees 

PRu=[ 1.05 2.07 4.78 4.48 
0.79 2.37 5.25 5.25 
0.65 2.57 5.47 5.47 
0.56 2.70 5.59 5.59 
0.51 2.80 5.66 5.66 
0.49 2.88 5.71 5.71 
0.49 2.95 5.75 5.75]'; 

%Polarization Ratio in dB 

Pr_Interpol_in_dB = interp2 (Spd_Wind_Tb, Beam_Inc_Tb, PRu, V, theta,... 

IntModeType); 
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B. CMOD-IFR2 IFREMER MODEL CODE FROM CCRS 


g. 

0 



g. 

o 


c 


X 


function sigmaO = ifremer(V,phi,theta); 

calculate radar cross-sections for cmod3 model wind retrieval 


algorithm according to IFREMER off-line analysis 
inputs 

V - the wind speed in m/s 

phi - the relative radar look-wind in degrees... this angle 
zero whe nthe wind is blowing towards the radar 
theta - the incidence angle in degrees 
output 


simgaO - radar cross-section (linear 


= [-2.437596 -1.567030 
-0.039714 0.410275 

-0.026785 0.130832 

0.192522 -0.027294 

0.003595 0.073938 

0.012663 0.044884 

-0.007864 -0.001635] 

= (theta - 36) ./19; 


0.370824 
0.184862 
-0.466007 
0.032854 
0.12845 
0.004448 


units) 


is 


pO = 1; 

pi = x; 

p2 = (3.*x.^2 - 1)./2; 

p3 = (5.*x.^2 - 3).*x./2; 

alpha = c(l) + c(2).*pl + c(3).*p2 + c(4).*p3; 

beta = c(5) + c(6).*pl + c(7).*p2; 

delta = c(8) + c(9).*pl + c(10).*p2; 

bO = 10.^(alpha + beta.*sqrt(V - delta)); 
bl = c(ll) + c(12).*pl + c(13).*V; 

b2 = c(14) + c(15).*pl + c(16).*p2 + (c(17) + c(18).*pl + c (19) .*p2) . *V 

+ c(20).*V.-2; 

sigmaO = b0.*(l + bl.*cos(phi.*pi/180) + b2.*cos(2.*phi.*pi/180)) 


C. INTERNET SITES: 

CCRS/RADARSAT, http://www.ccrs.nrcan.gc.ca/ccrs/tekrd/radarsat/rsate.html/ 
ENVISAT/ASAR, http://envisat.estee.esa.nl/instruments/asar/ 

ALOS/PALSAR, http://alos.nasda.go.ip/index-e.html 
JPL/NASA, http ://southport. jpi.nasa.gov/dese/imagingradarv3.html 
PIRATA, http://www.brest.ird.fr/pirata/pirataus.html 
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